Dynamic polarizabilities and related properties of clock states of ytterbium atom 
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We carry out relativistic many-body calculations of the static and dynamic dipole polarizabilities 
of the ground 6s 2 So and the first excited 6s6p 3 Po states of Yb. With these polarizabilities, we 
compute several properties of Yb relevant to optical lattice clocks operating on the 6s 2 — 6s6p 3 Pq 
transition. We determine (i) the first four magic wavelengths of the laser field for which the frequency 
' of the clock transition is insensitive to the laser intensity. While the first magic wavelength is known, 

, we predict the second, the third and the forth magic wavelengths to be 551 nm, 465 nm, and 413 nm. 

(ii) We reevaluate the effect of black-body radiation on the frequency of the clock transition, the 
resulting clock shift at T = 300 K being —1.41(17) Hz. (hi) We compute long-range interatomic van 
der Waals coefficients (in a.u.) C 6 (6s 2 ^So+Gs 2 %) = 1909(160), C 6 (6s 2 ^o+GsGp^o) = 2709(338), 
andC 6 (6s6p 3 Po+6s6p 3 Po) = 3886(360). Finally, we determine the atom-wall interaction coefficients 
(in a.u.), C3(6s 2 ^So) = 3.34 and C"3(6s6p 3 Po) = 3.68. We also address and resolve a disagreement 
between previous calculations of the static polarizability of the ground state. 
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PACS numbers: 06.30.Ft, 32.10.Dk, 31.25.-v 



INTRODUCTION 



Ytterbium atom is employed in a number of projects aimed at studying fundamental problems of modern physics 
by means of atomic physics. A large parity- violating signal on the 6s 2 "VSo - 6s5d 3 Z?i 408-nm transition of Yb has 
been recently observed at Berkeley A search for the CP-violating (T-, P- odd) permanent electric dipole moment 
of Yb was initiated Q ■ A Bose- Einstein condensate in a dilute gas of the ground-state Yb atoms has been attained Q 
C/3 ■ and a number of experiments with ultracold Yb atoms has been carried out (see, e.g., Refs. 0J,i). 

Here we focus on the optical lattice clock applications of Yb. In the optical lattice clocks, neutral atoms are trapped 
in a standing wave of a laser light (optical lattice) operated at a certain "magic" wavelength Q . At this wavelength a 
differential light perturbation of the two clock levels vanishes exactly Using Yb atom in lattice clocks was proposed in 
Ref. Q and experimentally demonstrated in Boulder [EHE]- Apart from numerous technical applications, such atomic 
clocks may be used for searching for a variation of fundamental constants [ll]| , mapping out atom- wall interaction [12] , 
, and serve as a basis for quantum information processing (see, e.g., (l3l. Il4|). 

In the present paper, we evaluate dynamic electric-dipole polarizabilities of the clock states, 6s 2 So and 6s6p 3 Pq . 
\ Dynamic polarizabilities of real and imaginary frequencies are linked to a variety of important atomic and molecular 
■ properties. In particular, we compute magic wavelengths, the effect of black-body radiation on the frequency of the 
clock transition, atom- wall interaction constants C3, and, finally, long-range molecular van der Waals coefficients Cq 
for the clock states. We take advantage of a close link between these properties and improve our theoretical predictions 
by referencing to measured properties, such as the first magic wavelength [E HE] and the van der Waals coefficient for 
two ground-state atoms @. 

Ytterbium has been studied both experimentally (see, e.g., [E [E [E HE HE HE HE HE] ) and theoretically (see, e.g., 
• j-h \ [1, HE HE HE HI, HE HI, HE HE]- Trapping and cooling schemes were discussed and implemented, the frequency of 
the clock transition was measured to a high accuracy [l(|,the first magic wavelength [E HE] and the van der Waals 
coefficient for the ground state of Yb were also measured [5| ■ 

Theoretical studies are not free from a controversy. Calculated polarizability of the Yb ground state reported in 
Ref. [20I ] is in a disagreement with most of other ab initio or semi-empirical analyses. The source of this disagreement 
is explained below and the discrepancy is resolved. The resolution is closely related to the role of the excitations 
from the core which are discussed in detail. Corresponding statement can be formulated as a general theorem: 
Consider a second-order transition matrix element, involving summation over a complete set of intermediate states. 
Then, a contribution from a subspace spanned by degenerate states does not depend on mixing of these states. This 
statement has implications not only for computing polarizabilities, but also for the two-photon transition amplitudes, 
parity-violating amplitudes, etc. in Yb and other atoms. 
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II. METHOD OF CALCULATION 



In this section we start by recapitulating expressions for dynamic polarizability and its relations to magic wave- 
lengths, and atom-atom and atom-wall interaction coefficients. Next we discuss our computational relativistic many- 
body method. 

The ac Stark shift of an energy of a spherically-symmetric, J = 0, atomic state v in electro-magnetic wave of 
amplitude e and frequency u) is given by 

A£,(a/) = -a»(|) 2 , (1) 
where the electric-dipole dynamic polarizability reads (we use atomic units in which H = 1, |e| = 1, and m e = 1) 

Here d is an electric dipole operator and summation is over a complete set of many-body states. 

Magic laser frequency uj* is found from the condition that the ac Stark shifts of both clock states v and w are 
identical so that the frequency of the transition does not depend on the laser intensity, 

AE v {uj) - AE w (uj) = -Mw) " «-M) (f) = 0. (3) 

Simply, at the magic frequencies, a v (u>*) = a w (u>*). 

Clock frequency shift due to black-body radiation is expressed in terms of the difference of static polarizabilities of 
the two clock levels (see detailed discussion in Ref. plj|). 

We also require atomic polarizabilities evaluated at purely imaginary frequencies 

These are used in evaluating the atom-wall interaction constant C3 (2?| 

1 f°° 

C3 = -7— \ a(iu)du), (5) 



47T _ 

and the van der Waals coefficients Cq. For two separated atoms in states w and v, 



Cq(w + v) = — / a w (iuj)a v (iuj)duj . (6) 
n Jo 

The Cq for two atoms interacting in identical states is subsumed by the above formula. Eq.([5]) is derived under 
assumption that there are no downward dipole transitions from either w or v state, which holds in our case. 

It is clear that for computing the polarizabilities, we need to determine electric dipoles and energies, and to perform 
the summation over a complete set of atomic states. To this end, we use a combination of the configuration-interaction 
(CI) method and the many-body perturbation theory (MBPT) to construct an effective Hamiltonian for two valence 
electrons (the CI+MBPT method [H, H^]). Further, we employ the Dalgarno-Lewis method 30] to reduce the 
summation over a complete set of many-electron states to solving an inhomogeneous system of linear equations. 



A. CI+MBPT method 

The ground state configuration of Yb is 4/ 14 6s 2 and most of the low-energy excited states correspond to configu- 
rations with one of the 6s electrons being excited. All these states can be represented to a relatively high accuracy 
as states with two valence electrons above a closed-shell core. This is the starting point of our calculations: the com- 
putational (CI) model space is spanned by all possible excitations of the two valence electrons. Implications arising 
from the states with excitations from the outer 4/ core sub-shell will be discussed below. 

The effective CI+MBPT Hamiltonian for two valence electrons has the form 



H cS = h 1 (r 1 ) + h 1 (r 2 ) + h 2 (r 1 ,r 2 ), 



(7) 
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where hi is the single-electron part of the relativistic Hamiltonian 

Ze 2 

hi = cap + - l)m e c 2 - + V^ -2 + (8) 

r 

and hi is the two-electron part of the Hamiltonian 

e 2 

h 2 {r u r 2 ) = - r +£ 2 (r 1 ,r2). (9) 

Fi _ r 2| 

In these equations, a and $ are the conventional Dirac matrices, V N ~ 2 is the Dirac-Hartree-Fock (DHF) potential 
of the closed-shell atomic core (N — 2 = 68, Z = 70), and £ is the correlation operator. It represents terms in the 
Hamiltonian arising due to virtual excitations from atomic core (see Ref. [HI, [29| for details) . £ = corresponds to 
the standard CI method. £i is a single-electron operator. It represents a correlation interaction (core-polarization) 
of a particular valence electron with the atomic core. £2 is a two-electron operator. It represents screening of the 
Coulomb interaction between the two valence electrons by the core electrons. We calculate £ in the second order 
of the MBPT. We use a B-spline technique 31] to construct a complete set of single-electron orbitals. We use 40 
B-splines in a cavity of radius R = 40 as and calculate the eigenstates of the V N ~ 2 DHF Hamiltonian up to the 
maximum value of the angular momentum l max — 5. The same basis is used in computing £ and in constructing the 
two-electron states for the valence electrons. 30 out of 40 lowest-energy states for every I up to l ma x — 5 are used 
to calculate £ and 14 lowest states above the core are used for every I up to l max — 4 to construct the two-electron 
states. 

The two-electron valence states are found by solving the eigenvalue problem, 

H cS V v = E V V V , (10) 

using the standard CI techniques. Calculated and experimental energies of a few lowest-energy states of Yb are 
presented in Table |IJ One can see that the pure ab initio energies are already close to the experimental values. 
However, for improving the accuracy further, we re-scale the correlation operator £1 by replacing £1 in the effective 
Hamiltonian ([7]) in each partial wave s,Pi/2,P3/2, • ■ ■ by / a £i- The rescaling factors are f s = 0.8772, f p = 1.03, 
fd = 0.933, and // = 1. These values are chosen to fit the experimental spectrum of Yb. The result of the fitting is 
shown in Table |TJ Note that the fitting is not exact for all the levels because the number of levels is larger than the 
number of fitting parameters. 



TABLE I: Energy levels of Yb (cm l ). The 4/ 13 5d6s 2 (7/2, 5/2)° level lies outside of the computational model space. 

State J Exp. [32] CI+MBPT 



ab initio Rescaled £1 









0.0 


0.0 


0.0 


4/ 14 6s6p 







17288 


18246 


17289 






1 


17992 


18946 


17996 






2 


19710 


20688 


19759 


4/ 14 5d6s 


:i D 


1 


24489 


24922 


24489 






2 


24752 


25195 


24743 






3 


25271 


25765 


25276 


4/ i4 6s6p 


Ipo 


1 


25068 


26463 


25611 


4/ 13 5d6s 2 


(7/2,5/2)° 


1 


28857 






4/ i4 5d6s 


L B 


2 


27678 


28485 


27812 


4/ 14 6s7s 


:1 S 


1 


32695 


33262 


32668 


4/ 14 6s7s 


l 8 





34351 


34871 


34280 



Transition amplitudes are found with the random-phase approximation (RPA) [33l |34| ] 

El vw = (* v \\dz + 6V N - 2 \\* w ), (11) 

where 5V N ~ 2 is the correction to the core potential due to core polarization by an external electric field. We compile 
our representative dipole matrix elements in Table [TTJ The values were computed in the CI+MBPT approach. We 
also compare with a high-precision result [l6| derived from fitting molecular long-range potentials to photoassociation 
spectra taken with ultracold samples. The 16% theory-experiment disagreement for the 6s 2 Sq — 6s6p P° amplitude 
is due to mixing of the 6s6p 1 P 1 ° state with the core-excited state 4/ 13 5d6s 2 (7/2, 5/2)°, which lies outside the 
computational model space (see Section Hill ) 
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TABLE II: Reduced matrix elements (a.u.) of the electric dipole operator for transitions in Yb. Values were computed in the 
CI+MBPT approach. The theory-experiment disagreement for the 6s 2 1 So — 6s6p 1 P° is due to mixing of the 6s6p 1 P° state 
with the core-excited state 4/ 13 5d6s 2 (7/2, 5/2)°, which lies outside the computational model space (see Section ILTT1 ) 



Transition 


Calc. 


Exp. 


6s' L S -6s6p a P? 


0.587 




Qs 2 'So - 6s6p 


4.825 


4.148(2)° 


6s6p 3 P ° -5d6s 3 L>i 


2.911 




6s6p 3 P ° -6s7s 1 S'o 


1.952 





Ref. mi 



B. Dalgarno-Lewis and RPA methods 



Computing polarizability requires summing over a complete set of many-body states. The CI+MBPT subspace of 
two valence electrons covers only a part of this set. The other class of intermediate state are core excitations, where 
the valence electrons serve as spectators. Notice, that even in the independent particle picture, the state of the valence 
electrons affects possible core excitations: core electrons can not be promoted into excited orbitals occupied by the 
valence electrons due to the Pauli exclusion principle. Following this discussion, we divide the polarizability into three 
parts [131, 

a(w) = aw (a;) + a COIe (u) + a coro _ V aiM . 

The valence contribution, a va x(u>), is given by where v and k are two-electron CI states (i.e., k lie entirely in 
the model space). We use the Dalgarno-Lewis method [30] for summing over a complete set of two-electron states. In 
this method, a correction S^ v to the two-electron wave function of the state v is introduced and the contribution of 
the valence electrons to the polarizability is expressed as (here we consider static polarizability) 

Oval(0) = !<**r,||d||*t,>. (12) 

The correction S^ v is found by solving the system of linear inhomogeneous equations 

(H eB -E v )5* v = -(d s + 5V N - 2 )9 v . (13) 
In case of dynamic polarizabilities we employ the identity 

Ek — E v _ 1 1 

{E k - E v f - cj 2 ~ (E k - E v ) - lo (E k -E v )+u' 

i.e., the equations (fT3|) need to be solved twice at two different energies E = E v ± lo. For polarizabilities of purely 
imaginary argument, lo — > iui in the above equations. 

The core polarizability is found using the RPA method [35] (linearized response theory). This method yields required 
eigen-energies of many-body particle-hole excitations from the core and relevant electric-dipole transition amplitudes. 
The dynamic core polarizability reads 

a COIC (iuj) = V - — . (14) 

uv>0 +^ 

Here the summation is over particle-hole excitations from the ground state of the atomic core; lo^ are excitation 
energies and / M are the correspon ding electric-dipole oscillator strengths. Technically, compared to the original, 
differential-equation method of Ref. [35| , we use a B-spline basis set expansion. An important property of the RPA core 
polarizability is that it satisfies (non-relativistically) the Thomas-Reiche-Kuhn sum rule, lim^—^ lo 2 a core (iLo) = N — 2. 



This property is especially important in evaluating the C3 and Cq coefficients for heavy atoms [27 1. 

Finally, the core- valence counter term, a coro _ va i(aj), is small. At the DHF level, we find that the relevant contribu- 
tions to static polarizabilities are below 1% and we neglect it at our present level of accuracy. 



III. THE ROLE OF THE 4/ 14 6s6p l P? - 4/ 13 5d6s 2 (7/2, 5/2)? MIXING 



There is a number of relatively low-lying states in the discrete spectrum of Yb which have an excitation from the 
outer 4/ sub-shell of the core. These states certainly cannot be considered as states with two valence electrons above 
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closed shells. One of such states is 4/ 13 5d6s 2 (7/2, 5/2)° which lies close to the 4/ 14 6s6p l P? state (see Table HJ. 
Notice that both states have the same total angular momentum and parity. Therefore the Coulomb interaction 
strongly mixes the two states. This mixing complicates the calculations. Our effective Hamiltonian ([7|) does take into 
account excitations from the core via second-order correlation operator £. Yet, the second-order approximation is 
not adequate in this particular case, as the two states are separated by mere 0.0173 a.u. (3789 cm -1 ). 

The calculation of the ground state polarizability is affected by the 4/ 14 6s6p 1 P° — 4/ 13 5d6s 2 6(7/2, 5/2)° mixing. 
The 6s 2 x So —> 6s6p 1 P° electric dipole transition contributes about 90% to the polarizability of the ground state. 
Corresponding theoretical CI+MBPT transition amplitude, 4.825 a.u., deviates significantly from the experimental 
value 4.148 a.u. As it was pointed out in Ref. [19] the reason for this disagreement is the above mentioned 

mixing. 

Calculation of the core polarizability displays a similar dependence on the mixing. Our RPA result for a core (0) = 
6.39 a.u. At the same time, semi-empirical result [24| based on the life-time measurements of low-lying core-excited 
states in Yb suggests a much higher value of 24(4) a.u. Again, the reason for this large difference lies in the strong 
mixing of the 4/ 13 5g?6s 2 (7/2, 5/2)° and 4/ 14 6s6p 1 P° states. Indeed, in the RPA core polarizability calculations, we 
consider a different system, Yb ++ , which is unaware of the presence of the two valence electrons. 

Since the 6s 2 4 So — ► 6s6p 1 P° transition contributes about 90% to the polarizability of the ground state, the 
authors of Ref. [2(| suggested to replace the CI+MBPT dipole matrix element by its experimental value. This was 
done to improve the accuracy of computing the polarizability. Below we argue that such a substitution cannot be 
justified unless a similar correction is done for the 4/7/2 — ► 5ds/2 transition amplitude in the calculation of the core 
polarizability. 

Let us use the following notation for the pure (unmixed) states 

$! = |4/ 14 6s6p X P?) , 

$2 = |4/ 13 5d6s 2 (7/2,5/2)°). 

The mixed states can be written as 

= $ x cos 9!) — $2 sm<f>, 
^2 = $1 sin (f> + $2 cos <f>, 

where <f> is mixing angle. 

Let also use ^0 = ®o for the ground state wave function. 

If we neglect a small energy difference between the states ^1 and \&2 (the relevant energy denominators entering 
a(0) differ by about 10%) then their total contribution to the polarizability of the ground state is proportional to 

(* ||d||*l) 2 + (* ||dp 2 ) 2 = 

($ ||d||$i) 2 + ($ ||d||$ 2 ) 2 (15) 

and does not depend on the mixing angle <\>. This means that either bases {<I>o, $2} or {^0, ^2} can be used 
in the calculations. The ab initio calculations within the CI+MBPT framework in which Yb is considered as a two- 
valence-electrons atom uses the non-mixed basis while the analysis based on experimental data naturally corresponds 
to the use of the mixed basis. Both approaches produce close results as will be demonstrated in Section |V] By 
contrast, the substitution of the experimental value for the (6s 2 1 So|D|6s6p 1 P°) transition amplitude is equivalent 
to using the mixed non-orthogonal basis {^Oj ^1, ^l- This seems to introduce a large error into the calculation of 
the polarizability. Indeed, the value reported in Ref. [20, HH is 111.3 a.u. which is about 30% smaller than the results 
of most ab initio calculations and semi-empirical analyses, compiled in Table Mil For example, the analysis [24| of 
the experimental data yielded aig o (0) = 136.4(4.0) a.u. 

The result of Ref. [2(| Hl[ is brought into an agreement with that of Ref. [24] by consistently correcting the core 
polarizability for the mixing. A DHF contribution from the 4/ subshell to core polarizability is 1.9 a.u. If we replace 
this value by the value based on experimental data (i.e., account for the mixing), 24(4) then the result of Ref. [2p( ] 
increases to mg o (0) = 133 a.u. consistent with the value of 136.4(4.0) a.u. from Ref. [24j . 

It is useful to find out what happens if the difference in energies of the two mixed states is not neglected. We 
will use a short-hand notation AE 1 = E(Vi) - E($ ), AE 2 = E($ 2 ) ~ E(V Q ), and S = 1/AE 2 - 1/AE^ Since 
the state * 2 = 4/ 13 5d6s 2 6(7/2, 5/2)° is above the state *i = 4/ 14 6s6p l P° the parameter 8 is negative. Further, 
A01 = ($ ||d||$i) and A 02 = ($ ||d||$ 2 >- 

A partial contribution of the two states into the polarizability reads 



5a = — — (Aq! + A\ 2 ) + S (A 01 sin (j> + A 02 cos<^) 2 . 



(16) 
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We see that when the core-excited states lie above the states with which they are mixed, the correction to the 
polarizability due to this mixing is negative. This is the case for both polarizabilities of Yb considered in this work. 
We would keep this in mind while analyzing the accuracy of the calculations in Section [V] 



IV. MIXING OF DEGENERATE STATES 



In the previous section we considered a particular case of mixing of the 4f 14 6s6p 1 P 1 ° and 4/ 13 5<i6s 2 6(7/2, 5/2)° 
states of Yb and demonstrated that this mixing can be ignored in the ab initio calculations of the ground-state 
polarizability of Yb. Notice that similar problems arise in computing the polarizability of the 3 Pq state and in many 
other cases not limited to calculations of the polarizabilities of Yb. Therefore, it may be useful to treat a general 
case and formulate the following theorem: Consider a second-order transition amplitude, involving summation over 
a complete set of intermediate states. Then, a contribution from a subspace spanned by degenerate states does not 
depend on mixing of these states. 

A generic contribution to a second-order transition amplitude can be written as (here Eq is some energy off-set) 

Mf = y {F f \P\F m )(F m \Q\F t ) 

Here Fi and Ff are some initial and final atomic states, and P and Q are some operators. Let's assume that the 
summation over the complete set of intermediate states F m subsumes a summation over a group of degenerate states 
T> and rewrite the corresponding partial contribution as 

AA *f ^ E ( F f\P\ F ™){Fm\Q\Fi)- (18) 

Here the summation ranges only over degenerate states and AE is their common energy denominator. 

As in the previous section, we introduce the pure non-mixed basis states and rewrite the real (mixed) states F m 
as an expansion 



F m = ^2c mk <P k . (19) 



Both sets $ and F can be made orthonormal, with the matrix C — {c m k} being a unitary matrix, i.e., CC^ = I. 
Substitution of GUI into dTHJ leads to 

AE- 

kl 



kl m 

Here we used the ortho- normality condition, Cfe m c^ = 5 m i. We see that the final expression does not depend on 
the mixing coefficients c. 

We stress two conditions under which this theorem is valid. It is assumed that the unmixed states possess the 
same energy and that there is no mixing with the states outside the degenerate subspace. These two conditions are 
probably never fulfilled exactly. However, there is a large number of cases when they are fulfilled approximately. One 
of the criteria is that the average energy spread in the quasi-degenerate subspace \SEx>\ <C \E n — E-p\, where Ex> is 
the average energy in the subspace. 

Why is this practically important? In Yb, to rigorously account for the mixing, one has to extend the two- 
valence-electron model space to include computationally-expensive core-excited states. The theorem claims that 
while computing polarizabilities, a much smaller two- valence-electron space would suffice. In other words, at the CI 
stage, the excitations from the 4/ subshell can be safely ignored (even) if there are nearby two-electron states of the 
same symmetry. (Of course, virtual core excitations are included in the self-energy MBPT operator £.) Excitations 
from the core will have to be included into calculation of the polarizability of the core. This corresponds to working 
with the non-mixed basis 

The independence on mixing explains why pure two- valence ab initio calculations of polarizabilities give good results 
for Yb despite its complicated structure. This should be true not only for polarizabilities but also for two-photon 
transition amplitudes, parity non-conservation, etc., and not only for Yb but for some other atoms as well. 
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V. RESULTS AND DISCUSSION 
A. Static polarizabilities 

Our computed static polarizabilities of the clock states are presented in Table ITTll There we tabulate values obtained 
in various approximations and we also compile results from the literature. Below we analyze our results and estimate 
theoretical uncertainties. 

The first line of Table IIIII lists the results of pure ab initio calculations, next line gives the results obtained when 
energies are fitted as explained in Section [ill Then, to gauge the accuracy of the calculations, we introduce various 
corrections to the polarizability of the ground state based on available experimental data. In the third line of the Table 
we correct to reproduce the experimental value for the magic frequency 0, [l5| . This correction will be explained 
in the following section. The fourth value (135.2) is obtained by replacing the dominant ab initio contribution by 
semi-empirical values. This corresponds to changing from the unmixed $i, $2 basis to the mixed ^i, ^2 basis (see 
Section [Till for details). Specifically, we use the value 4.148 a.u. [l6| for the (6s 2 1 So||d||6s6p 1 Pi) transition amplitude 
and 5a = 24(4) [24| for the contribution of the excitations from the 4/ core subshell. The resulting value is smaller 
than both ab initio results (with and without fitting) consistent with our observation made at the end of Section IIIII 
Finally, we rescale a(0) of the ground state to fit the experimental value for the van der Waals coefficient Cq |5j. Since 
Cq is simply obtained by integrating a(iuj) 2 , Eq. ([5]), and both the integral and the static polarizability are dominated 
by the principal transition, the uncertainties of a(0) and of Cq are correlated (see a detailed discussion in Ref. [36j). 
The final value for the ground state and its uncertainty, 141(6), are chosen to cover the spread of our numbers in the 
Table. 

Now we turn to the polarizability of the 3 P°> state. The central point of our final value is chosen to coincide with 
the fitted value. Choosing the lower calculated value as the central point is justified by the fact that there is a 
correction due to excitations for 4/ and this correction is negative (see the end of Section IIIII for details) . The value 
of this correction must be smaller than that for the ground state because corresponding states with excitations from 
4/ are higher in the spectrum and closer to the neighboring two-electron states of the same symmetry. Notice that 
the polarizabilities of both clock states end up having a comparable fractional accuracy. Qualitatively, this may be 
explained that while the convergence with respect to increasing basis is slower for the 3 Pg state, calculations for the 
1 So state are affected stronger by the core-excited states. 

Results of other calculations and semi-empirical analyses are presented in Table HTT1 for comparison. For the ground 
state note a good agreement with the value of Ref. [Hj], a(0) = 136.4(4.0), derived using only experimental data. The 
value a(0) = 143 recommended by the same authors is also consistent with our result. 

The result for the ground-state polarizability is in a disagreement with previous calculations [H, HE Hl|. As 
discussed in Section [TTT1 the disagreement for 1 So comes from the inconsistent use of experimental matrix element for 
the principal transition in these works. As for the 3 Pg state, the present result and that of Ref. [2l| differ by about 2 
standard deviations and in experimental work this level of agreement between two values would be acceptable. 

Finally, let us emphasize that evaluating accuracy of theoretical calculations is a non-trivial exercise and due to 
the lack of high-accuracy experimental data for the 3 Pq state, the uncertainty in Ref. [2l[ was estimated as a half of 
the difference between pure ab initio result and the result obtained with the fitting of the energy. This may be an 
unreliable estimate. Suppose we carry out a similar analysis based on the present calculations for the ground state. 
Fitting energy in Si moves our ab initio result, 138.9, to 145.7, so using prescription of Ref. [2l[ would yield 145.7(3.4) 
which is displaced compared to and is as twice as accurate as our final value of 141(6). 

B. Magic frequencies 

Magic frequencies are frequencies of the laser field at which the ac Stark shifts of both clock levels are the same so 
that the frequency of the clock transition is insensitive to the laser intensity. Magic frequencies, uj*, are found from 
the condition 

m s >*) = a» P >*). (21) 

First five calculated magic frequencies for Yb are presented in Table IIVI The first frequency was measured to a 
high precision 0, [l5j]. Our ab initio calculations reproduce it with the 1.3% accuracy. The accuracy for the other 
frequencies is likely to be worse. This is because the energy denominators for the dynamic polarizabilities ^ are 
shifted by the frequency of the laser field and the assumption that the difference in energies of the 4/ 14 6s6p 1 P 1 ° 
and4/ 13 5c?6s 2 6(7/2, 5/2)° states can be neglected becomes less accurate (see Section UTTI for details). 

To improve the accuracy of predicting the magic frequencies and corresponding polarizabilities we modify formula 
@ in two ways. First, we use experimental data for the dominant terms: (6s 2 1 So||c?||6s6p 1 Pi)=4.148 [16j and 
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TABLE III: Static electric-dipole polarizability a(0) (in a.u.) of the 6s 2 1 So and 6s6p 3 Po clock states of Yb atom in various 
approximations and comparison with literature values. 





a(0), 3 Pg 


Comment 


this work 






138.9 


315.9 


ab initio with no fitting 


145.7 


301.6 


energies are fitted by rescaling Ei 


139.1 




with a correction to a( 1 So) to fit the magic frequency u — 0.0600 a.u. 


135.2 




using experimental data for dominant contributions 


138.8 




rescaled to fit experimental data[5] for Cq — 1932(30) a.u. 


141(6) 


302(14) 


final 


other 






131.6 




Wang et al [3JJ; Density functional theory, 1995 


145.3 




Wang and Dolg [38]; CCSD(T), 1998 


1/11 *7 

141. f 




Miller |oy]; IJensity lunctional tneory, zuUz 


118(45) 


252(25) 


Porsev et al [19]; CI+MBPT, 1999 


111.3(5) 


266(15) 


Porsev and Derevianko [20,21]; CI+MBPT, 2006 


154.7 




Buchachenko et al [40J; CCSD(T), 2006 


136.4(4.0) 




Zhang and Dalgarno [24]; based on experimental data, 2007 


157.3 




Chu et al [23]; Density functional theory, 2007 


143 




recommended value by Zhang and Dalgarno [24], 2007 


144.59 




Sahoo and Das [261; relativistic coupled-cluster theory, 2008 



(6s 2 1 So||(i||4/ 13 5(i6s 2 (7/2, 3/2)°)=2 [4l| . Second, we introduce a correction to simulate the effect of other excitations 
from the 4/ core shell, 

Aa = — ( + -LS) , (22) 

where E = 37414.59 cm -1 is the energy of the nearby state with the excitation from the core and A is a fitting 
parameter. Its value is chosen to fit experimental magic frequency. The rest of contributions are taken from the ab 
initio calculations. 

The above empirical correction has been introduced for the ground state only. The polarizability of the upper 
clock state was computed within the CI+MBPT method for the valence contribution and the DHF method for core 
polarizability. 

The values of magic frequencies and corresponding polarizabilities calculated with this empirical correction are 
listed in the third and forth columns of Table IIVI and the ac polarizabilities are shown in FigfTJ We see that the 
results for the first four magic frequency change little, while those for the fifth magic frequency vary significantly. 
Still, while difficult to predict, the fifth magic frequency may be of a practical interest for designing better clocks: by 
contrast to the first four frequencies, the ac polarizability here may be negative (that was the ab initio prediction, but 
the empirical correction flipped the sign of the polarizability). Because of that, the atoms are trapped at the minima 
of the laser intensity, which reduces perturbations of the clock frequency [42l |. 

Finally, using the modified formula at to = provides another accuracy test for our predicted static polarizability 
of the ground state (Section [V Ap . The resulting value, 139.1 a.u. was listed in Table [TTTl and used in estimating the 
accuracy of the calculations. 



C. Black-body radiation shift 

The effect of black-body radiation (BBR) on the frequency of the Yb clock transition is expressed in terms of the 
static polarizabilities. The relevant clock shift reads (see detailed derivations in Ref. [2l|) 

Svbbr « -^Kfs^) 3 T 4 (a 3Po (0) - m So (0)) , (23) 

where T is the temperature and «hfs is the fine-structure constant. With our polarizabilities we find at T = 300 
K, Svbbr = —1.41(17) Hz which is in good agreement with the result, Svbbr = —1.34(13) Hz, of Ref. [2l| despite 
our significant differences in values of individual polarizabilities. The agreement is fortuitous, as both polarizabilities 
of Ref. [2l| are smaller, most of the shift canceling out while taking the difference in Eq. (|2"3"|) . Note, however, 
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FIG. 1: (Color online) Dynamic polarizabilities of the two clock levels. "Magic" conditions occur when the two polarizabilities 
intersect off resonance. These are marked by small circles on the plot. 



TABLE IV: Magic wavelengths for the 6s 2 ^So — > 6s6p s Po transition in Yb and values of dynamic polarizabilities at corre- 
sponding frequencies (a.u.). 



A,nm 


ab initio 


a(u>),a.u. 


A,nm 


Corrected" 


a(oj),a.u. 


Experiment 
A,nm 


749.0 




193 


759.37 




186 


759.355 6 


552.2 




46 


551.5 




60 




459.3 




478 


465.4 




382 




413.2 




1246 


413.25 




817 




359.7 




-742 


402. 55 c 




365 c 





n aig (u;) is corrected to fit experimental magic frequency 

6 Ref. HI 

Unreliable 



that a somewhat smaller uncertainty reported in |2lj is mostly due to overestimating the theoretical accuracy for 
polarizabilities, especially for the ground state. See Sections IIIII and [V Al for details. 

It is worth emphasizing, that at the time of writing, the BBR correction is the leading uncertainty in the error 
budget of the lattice clocks [ToJ . The present calculations demonstrate a difficulty of atomic theory in determining 
this important correction accurately. We hope that our analysis would motivate a further experimental work, for 
example, on cryogenically cooled clock chambers to reduce (and ultimately to accurately determine) the BBR shift. 



D. Atom- wall interaction coefficients Cz and van der Waals coefficients C(, 

We calculate the atom-wall interaction constants C3 and the van der Waals constant C§ for both clock states using 
the formulae (O and © . The results are presented in Table [V] together with the results of other calculations and 
measurements for the Cg coefficient for the ground state. 

For two Yb atoms in their ground states we find Cq^Sq + 1 S ) = 2041 a.u. (CI+MBPT with rescaled £ to fit the 
energies). However, from our calculations of the polarizabilities we know that the contributions from the excitations 
from the 4/ core state are likely to make the results smaller. Therefore, we rescaled the central point for C^Sq + 1 SVj) 
at twice the rate as for the polarizability of the ground state (i.e., reduced it by 6.5%). This moves the the central 
value to 1909 a.u. Since the Cq is obtained as a quadrature of the square of polarizability, we assign a fractional 
uncertainty to Cq as twice as large as that for the static polarizabilities. Note a good agreement of the ground-state 
Cq with recent experimental determinations || and with calculations (24j. The results in the Table for Cq[}Sq + 3 Pq) 
and Cq( 3 Pq + 3 Po) were obtained in CI+MBPT with rescaled £ to fit the energies. 

Including core polarizability into computation of Cq is important. We find that neglecting a COIC (iuj) reduces the 
final result for Cq by about 10%. This is consistent with the earlier observation (2?| for heavy alkali-metal atoms. 

Atom-wall interaction coefficients C3 for the clock states were obtained in CI+MBPT with rescaled £ to fit the 
energies; these are also listed in the Table [V] Notice, that the difference in C3 coefficients for the clock states was 



10 



computed by us earlier in Ref.|l2|]. While forming the difference the contribution of core polarizability cancels out. 
For individual states, this contribution is substantial. For example, for the ground state result is increased from 2.09 
to 3.34 a.u. due to a COTe (iuj). As a large fraction of C3 value is accumulated due to the core excitations, assigning 
error to C3 is difficult For example, to estimate error in C3 for simpler alkali-metal atoms in Ref. [27| , results of two 
methods of computing C3, (i) via the integral of a(iu)) and (ii) by computing expectation value of some two-body 
dipole-dipole interaction were compared. A highly-technical evaluation of the expectation value of two-body operator 
is beyond the scope of the present paper and we do not assign uncertainties to C3. We expect that the errors in C3 
are unlikely to exceed 50%. 



TABLE V: Atom-wall interaction coefficients C3 and van der Waals coefficients Ce (a.u.) for the 6s 2 ^Sb and 6s6p 3 Po states of 
Yb. 





Value 






Comment 






Atom-wall interaction 






c 3 


'So 


3.34 


Calc. 


this work 


c 3 


3 Pg 


3.68 


Calc. 


this work 






Atom-atom interaction 






Ce 


So + 1 So 


1909(160) 


Calc. 


this work 






2300(250) 


Exp. 


Enomoto et al [171 






1932(35) 


Exp. 


Kitagawa et al [5] 






2400 - 2800 


Calc. 


Buchachenko et al [40J 






2292 


Calc. 


Chu et al [23] 






2062 


Calc. 


Zhang and Dalgarno [24] 


C e 


' A PS+ 3 PS 


3886(360) 


Calc. 


this work 


C e 


1 s + 3 ps 


2709(338) 


Calc. 


this work 



VI. CONCLUSION 

We calculated polarizabilities of the clock states of Yb with the accuracy of about 5%. The substantial disagreement 
with previous calculations for the ground state polarizability was explained and resolved. We also computed the 
first four magic frequencies of the lattice laser field, the effect of black-body radiation on the frequency of the 
clock transition, the C3 atom-wall interaction constants and the Cq van der Waals coefficients for both clock states. 
Polarizabilities, the first magic frequency and the Cq coefficient for the ground state are in a good agreement with 
the most accurate calculations and measurements. The presented data may be of interest for designing better clocks, 
applications of the clocks in studying atom-wall interaction and quantum information processing, and quantifying 
molecular potentials for ultracold collision studies. 
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